A New Mechanism for Bubble Nucleation: Classical Transitions 
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Given a scalar field with metastable minima, bubbles nucleate quantum mechanically. When 
bubbles collide, energy stored in the bubble walls is converted into kinetic energy of the field. This 
kinetic energy can facilitate the classical nucleation of new bubbles in minima that lie below those of 
the "parent" bubbles. This process is efficient and classical, and changes the dynamics and statistics 
of bubble formation in models with multiple vacua, relative to that derived from quantum tunneling. 



Consider a potential V((f>) with many metastable min- 
ima with positive vacuum energy. A typical region in 
the universe undergoes de Sitter expansion, with Hub- 
ble parameter H ~ y/V(<pi)/M p , where M p is the usual 
reduced Planck mass and <f>\ labels the vacuum. Small 
regions may tunnel to another minimum, 

02, V(<j> 2 ) < 

V{4>i) [I] , forming an expanding "bubble". Some of these 
bubbles inevitably collide. In this Letter we show that 
these collisions can give the field sufficient energy to climb 
over potential barrier (s) and form new bubbles at min- 
ima where V(<p) is less than that of cither of the original 
bubbles. Bubble collisions can thus yield new bubble 
universes - the progeny become the protagonists. In a 
universe with many metastable vacua this new classical 
mechanism populates the de Sitter sea with additional 
bubbles, beyond those provided by tunneling. 

Such classical transitions should not be surprising. 
Hawking, Moss and Stewart |2J numerically analyzed col- 
lisions in a potential with two minima. When cosmolog- 
ical bubbles of the lower vacuum collide, they can form 
a new bubble of the higher, metastable vacuum. This 
new bubble is surrounded by regions in the lower vac- 
uum state and thus quickly collapses. Similar features 
are seen in numerical simulations by |3J0]. The key new 
result here is that, given a third, lower energy, minimum, 
the field can "slosh" into that vacuum after a collision 
between bubbles in higher energy minima. This yields a 
new bubble which is stable and expanding, since its inte- 
rior energy density is lower than that of the surrounding 
region. 

Physically, this process is easy to understand, even 
though the detailed dynamics require numerical simula- 
tions. A bubble collects energy as it expands into a region 
of higher potential energy, storing it as the gradient en- 
ergy of its walls. When two bubbles collide, this energy 
must go somewhere. The bubble wall may lose coherence 
and radiate the energy into the space in which the two 
bubbles have merged [5]. Alternatively, the energy in the 
walls can be deposited into kinetic energy (in field space) , 
allowing the field to climb potential barriers and thus 
slosh into adjacent minima. If one of these minima has a 
lower potential energy than the parent bubbles, we find 



that this process is efficient, provided the parent bubbles 
have had a (small) amount of time to expand since they 
were nucleated, allowing them to store sufficient energy 
in their walls. Clearly, the transition condition depends 
on the size of the barriers in the model, relative to the 
energy stored in the walls. After nucleation, the collision 
energy is deposited into the new bubble wall, hence the 
new bubble is formed with non-stationary walls. 

Transition Condition: Assume a setup illustrated in Fig. 
[T] with three local minima 4>ii4>2i<t>3- Two bubbles of </> 2 
are quantum mechanically nucleated within the false vac- 
uum of 0i. Subsequently, the two </>2-bubblcs collide, and 
the question is: under what condition would the collision 
trigger the classical nucleation of a third bubble of 03? 
For an analytic estimate, we will work in the flat space 
limit. The hyperbolic foliation is convenient for studying 
coUisions 0[71[H]: 

ds 2 = -dt 2 + dx 2 + t 2 (diP 2 + swh^dO 2 ). (1) 

The two 02-bubbles are nucleated at t = 0, x = ±6 (and 
?/' = 0). The pre-collision evolution is well-described by 
the solution of pQ. We make the thin- wall approximation, 
where each bubble wall moves according to Rq — (x ± 
b) 2 — t 2 , with R being the initial bubble radius. The 
wall thickness at a given t follows Ax = 7~ 1 Axo, where 
Axq is the initial thickness, and 7 is the Lorentz factor: 
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Assuming b Rq, the collision occurs at x = where the 
Lorentz boost is 7 = b/Ro (hereafter 7 takes this value). 
The initial wall thickness can be estimated from the 
Compton wavelength i.e. Ax ~ A(pi 2 / y/ AV^ 2 , where 
A0i2 = <fii — 4>2 and AV^ 2 is the barrier in between. 

The post-collision evolution follows from □(/> = d<f,V , 
which in the coordinates of Eqn. (JlJ is 

d 2 (j) 2d(f>_d 2 (f> dV 
1W + tlk ~ dx? ~ 1m> 



(3) 



Let us focus on the time evolution of <f> at x = (and 
implicitly ip — 0; see Fig. [2]). Pre-collision, <f> is stuck 
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at the false vacuum of 0i, where dV/d(f> = 0, and both 
spatial and time derivatives of vanish. As the collision 
begins, is still at </>i, hence the potential force dV/d<f) 
remains zero, but the presence of the two bubble walls 
means d 2 (f>/dx 2 is non-zero and gives rise to a kick to- 
wards 02 (i-e. negative in the convention of Fig. [IJ: 

( 4 ) 

ftc 2 7 A0i 2 ' 1 ' 

where we have used the 1-bubble results mentioned ear- 
lier. This kick initiates a field motion from 0i tending 
towards 02- Whether this motion can be completed de- 
pends on whether the kick wins over the counteracting 
potential force (~ Ay i b 2 /A0 12 ) that inevitably develops. 
Absent dissipation, to which we will return below, Eqn. 
Q suggests the kick always wins since we expect 7 > 1 
in general. 

Once the field makes it over the first barrier between 
vacua 1 and 2, it would appear the field can also overcome 
the second barrier between vacua 2 and 3, if the second 
barrier is smaller than the first (absent dissipation). If 
not, Eqn. ^ and simple energy conservation suggests 
the required condition for overcoming the second barrier 
is 7 2 > AV^/AV^ (see Fig. [j] for definitions). Sum- 
marizing, the condition for a collision induced classical 
transition from vacuum 1 to vacuum 3 is 

7 2 >— max(l,^|y (5) 

1 1-13 V V ; 

Here, we have introduced a factor of 1/(1 — (3) to approx- 
imate the effect of dissipation, where (3 can be thought 
of as the fractional energy dissipation. There could be 
many sources of dissipation: the second term on the left 
of Eqn. ([3]) is one example; other examples include radi- 
ation into fields is coupled to (including itself), grav- 
itational waves (if there are significant deviations from 
spherical symmetry), and Hubble friction. The last item 
is simple to estimate: Hubble friction becomes impor- 
tant if A0i3 > H~ x d<j>/dt ~ H-^A^u/Axo, where 
A0i3 = 0i — 03. Since the initial bubble wall thickness 
Axo is generally a small fraction of H~ x , we expect Hub- 
ble friction to be unimportant unless A0 13 3> 7A0 i2 . 

If Eqn. ([5| is satisfied, the (model-dependent) dissipa- 
tion likely causes the field at x = to eventually settle 
in 03, that is, unless there are additional lower minima 
to the left, in which case further field excursions are pos- 
sible. It is interesting to note that Eqn. ([5| can always 
be satisfied if the collision is sufficiently relativistic, i.e. 
if enough time elapses between nucleation and collision. 
If j3 is not too close to unity, and if AV^ 2 ~ AF 13 , clas- 
sical transitions over potential barriers are generic even 
if the collision is barely relativistic. On the other hand, 
if Eqn. (|5| is not satisfied, the two bubbles merge and 
the energy stored in the walls dissipate via one or more 
of the dissipation mechanisms outlined above. 
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FIG. 1: Schematic diagram of a 3-minima potential. The 
false vacuum at 0i quantum mechanically nucleates bubbles 
into regions of 02 • Sufficiently energetic collisions of two 02- 
bubbles can classically nucleate a bubble of 03. The two bar- 
rier heights are defined with respect to V(cf>i). Superscript 'b' 
means barrier. 

We numerically test these ideas using a toy model: 

V(4>) = ^0 2 (0 2 - 2 ) 2 + eA0^(0 + O ) + a \4. (6) 

We assume e is small enough so that we have three clas- 
sically stable points, 0i 0o, 02 ~ and 03 — 0o- 
The barriers are parametrized by 0o an d A, while a sets 
the overall vacuum energy scale which we assume is much 
larger than the difference in energy between the minima. 

These minima are perturbatively non-degenerate, and 
we can use the analytic solution of [1] as an approximate 
solution to tunneling events. Two bubbles of = 02 
nucleate within a sea of highest metastable vacuum = 
0i, with an initial radius i?o = (v'32Ae0 2 1 /3)~ 1 . These 
bubbles have an approximate initial profile, 

1/2 

0(r) = 0o (l + 2e-^° ir - R ^ , (7) 

where r is the distance from center of the bubble. The 
approximate height of the barrier between the top two 
minima is 

A^-f, (8) 

which is also roughly the barrier between the middle and 
lower minima AV^ ss AV^. We explore the transition 
condition numerically in three spatial dimensions, using 
a modified version of LatticeEasy [9] and 1024 3 lat- 
tices. We take the overall vacuum energy to be GUT 
scale, cvA0o w (10 _6 m p i) 4 and e = 1/30. We include a 
homogeneous expanding background. 

We begin with a case where the transition does not 
occur. The two bubbles nucleate close together, at a 
separation 2.4i?o- I n a static universe, the bubbles would 



FIG. 2: Conformal diagram showing the values of <j> on the X — r plane (left) and the Y — t plane (right). Our conformal 
coordinates, (t, X, Y, Z), are chosen such that X is the axis through the centers of the bubbles and r is defined in the caption 
of Fig [3] The collision occurs at X = (analogous to x — in the text). Regions of (fix appear off-white, regions of appear 
grey and regions of <f>3 appear black. In this simulation, the centers of the two bubbles are ~ apart at nucleation. 

When these bubbles collide (at r w 17.5) the walls have a Lorentz factor 7 w 1.5. The slight ripples in the middle of each 
bubble are numerical artifacts. 



collide when r = 1.2Rq, but expansion of the background 
delays the collision until r = 1.3i?o and 7 ss 1.3. Figure[3] 
shows the time evolution of this scenario. The field does 
not have sufficient kinetic energy after the collision to 
climb over the potential barrier, and hence the bubbles 
merge into a single region with <f> = <f>2 w 0. 

Now consider two bubbles which nucleate at a sepa- 
ration of 3Rq. The bubble walls achieve 7 « 1.8 before 
collision, and then successfully traverse the potential bar- 
rier. Figure [2] shows a conformal diagram of the resulting 
field profile, showing the classical nucleation of a bubble 
with vacuum energy V{(j)^). We ran several simulations, 
varying the distance between nucleation points. Generi- 
cally, we find that the critical Lorentz factor at which a 
transition occurs is 7 ~ 1.4, or fj « 0.5 via Eqn. ([5]). 
There are several interesting features. First, classical 
transitions occur even when the bubbles are only mildly 
relativistic. Second, if transitions do not occur, most of 
the energy is then released as debris (radiation) after the 
collision. Third, the collision induced bubble has a very 
small initial size (of the order of wall thickness) and is 
nucleated with a non-stationary wall. The bubble is ho- 
mogeneous inside and initially non-spherically symmet- 
ric, though later expansion tends to make it more spher- 
ical. Fourth, our 3D simulations confirm the robustness 
of the hyperbolic symmetry: grossly unstable symmetry 
violating perturbations are not seen. 

Discussion: We have shown that, for a single scalar 
field model with several minima, bubble collisions gener- 
ically lead to the formation of new, lower vacuum bub- 
ble unless the orignal bubbles nucleate very close to one 
other. This mechanism creates a new set of possibilities 
for old problems. 

1. Are classical transitions the dominant bubble forma- 
tion process ? Consider again the setup in Fig. [I] Let the 
nucleation rate from 1 to 2 be T12, and that from 2 to 3 
be r 2 3. If r i2 is too high, ^-bubbles will form too close 
to each other, and hence their wall energies will be too 
small for them to classically transition to 03-bubbles dur- 
ing collisions. In this case the universe percolates rapidly 



into the 4>2 vacuum, which eventually quantum nucleates 
03-bubbles. For this not to happen, the (^-bubbles need 
to be nucleated at a typical separation, ~ r i2 1 ^ 4 , that is 
larger than the initial bubble radius Rq by a minimum 
Lorentz factor of 7 m i n . (7min ^ s defined by the right hand 
side of Eqn. [5J). Therefore, a necessary condition for 
classical transitions to be important is 

r i2 < (7min.i? )" 4 (9) 

Whether classical transitions dominate over quantum 
tunneling into <j)$ depends further on T23. The question 
involves complex measure issues: given a random point 
in the <p3 vacuum, what is the probability that it arose 
from quantum tunneling or classical collisions? We will 
address this question elsewhere. 

2. Implications for the eternally inflating stringy land- 
scape. Collisions, if they are sufficiently relativistic, pro- 
vide a new way of scanning an eternally inflating land- 
scape. A collision could allow classical transitions not 
just over one or two barriers but multiple barriers. It is 
thus necessary to revisit predictions for bubble counting 
measures based on quantum tunneling alone [TOl E] • We 
caution however that multidimensional potentials gencr- 
ically possess complicated intra- field couplings, and such 
couplings can change the dynamics of the collision. For 
example, consider our toy model with an additional cou- 
pling to a light field \- This light field can be excited 
in a collision, and carry away energy thus increasing j3, 
making the collision less elastic. 

3. Observational signatures? It is conceivable the bub- 
ble which is our own universe is formed via a classical 
transition of the type discussed here. Our simulations 
suggest the bubble is quite homogeneous but is initially 
highly anisotropic. Is the subsequent expansion sufficient 
to make it acceptably isotropic? Or do we need some pe- 
riod of slow-roll inflation to make it both isotropic and 
flat? Is the residual anisotropy observable? The observa- 
tional signatures are likely different from those considered 

by pen en. 
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FIG. 3: Time evolution of two bubbles whose centers are ~ .36H 1 apart at nucleation. When these bubbles collide (at r « 9.0) 
the walls have a Lorentz factor 7 ~ 1.3. We use a conformal time a(t)dr = y/\4>odt, where a(0) = 1 at the beginning of the 
simulation and dt is the usual proper time of an FRW spacetime. Notice the Lorentz contraction of the wall thickness as the 
bubble expands. 
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FIG. 4: Time evolution of two bubbles whose centers are ~ A5H 1 apart at nucleation. When these bubbles collide (at 
t ~ 17.5) the walls have a Lorentz factor 7 = 1.8. We use a conformal time a{t)dr = V^cfi^dt, where a(0) = 1 at the beginning 
of the simulation and dt is the usual proper time of an FRW spacetime. 



4- Implications for the small cosmological constant 
problem. Abbott [13] (see also [14]) proposed a model 
for relaxing the cosmological constant using a step-wise 
potential via a series of tunneling events. Collisions in- 
troduce new, and perhaps faster, excursions through the 
multiple descending vacua. Tapping the collision energy 
might even help alleviate the well known empty universe 
problem [15] . 

To summarize, we presented a new mechanism of bub- 
ble nucleation where a potential barrier (or even multi- 
ple barriers) can be transitioned classically via energy re- 
leased during collisions of two or more bubbles. This new 
mechanism possesses a rich phenomenology which we will 
explore, along with improved analytic descriptions of the 
bubble formation criteria, in forthcoming publications. 
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